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■ , Abstract 

^ • In single-bubble sonoluminescence, a bubble trapped by a sound wave in a flask of 



liquid is forced to expand and contract; exactly once per cycle, the bubble emits a very 
^ ! sharp (< bOps) pulse of visible light. This is a robust phenomenon observable to the naked 
■ eye, yet the mechanism whereby the light is produced is not well understood. One model 
that has been proposed is that the light is " vacuum radiation" generated by the coupling of 
the electromagnetic fields to the surface of the bubble. In this paper, we simulate vacuum 
radiation by solving Maxwell's equations with an additional term that couples the field to 
the bubble's motion. We show that, in the static case originally considered by Casimir, we 
reproduce Casimir's result. In a simple purely time-dependent example, we find that an 
instability occurs and the pulse of radiation grows exponentially. In the more realistic case 
of spherically-symmetric bubble motion, we again find exponential growth in the context 
of a small-radius approximation. 



I. Introduction 

Single-bubble sonoluminescence [1] is a mysterious phenomenon. A small bubble of 
gas, usually air, is trapped at the center of a flask of liquid, usually water, by the application 
of an intense acoustic field. The frequency of the field is typically 25 or 30 kHz, and 
once per cycle, driven by the sound field, the bubble undergoes expansion and then rapid 
contraction. If the parameters are right [2] (here "parameters" refers to such things as 
the intensity of the sound field, the concentration of the gas, the chemical composition of 
the gas, the temperature of the water) the bubble will emit a very narrow pulse of light 
during the contraction phase of the cycle. The sound wave, with a time scale of tens 
of microseconds, produces a contraction of the bubble measured in tens of nanoseconds, 
which in turn somehow generates a pulse of visible light whose duration has recently been 
measured to be tens of picoseconds [3]. Furthermore, even though the motion of the bubble 
is quite violent, if the parameters are right it can be remarkably stable, repeating itself 
over millions of cycles, with the flash of light appearing at the same point in the cycle each 
time [4]. 

Other aspects of the phenomenology of sonoluminescence are also worthy of note. 
For example, the spectrum of emitted light is only partially known, because the water 
absorbs all wavelengths shorter than about 180 nm [5]. The part that is observed looks 
like the tail of a rising distribution, and attempts to flt it to a thermal spectrum have led 
to the speculation that the emitting region is very hot, certainly in excess of 25,000 degrees 
Kelvin, and perhaps even as high as a million degrees, at which point nuclear fusion might 
be expected to be significant [6] . 

Another peculiarity is the fact that, whereas air- filled bubbles work well as a vehicle for 
sonoluminescence, bubbles filled with either oxygen or nitrogen, or indeed with a suitable 
mixture of these two gases, do not [7]. The small noble-gas component of air is essential for 
significant sonoluminescence to take place. This agrees with the suggestion that, during 
the first second or so of the bubble's oscillation, the oxygen and nitrogen are ionized and 
absorbed by the water, leaving a rarefied bubble filled with noble gas [8]. Experiments done 
with bubbles filled with various noble gases confirm that they produce sonoluminescence 
efficiently. 

There is much more to the phenomenology of sonoluminescence, as described in a 
number of recent reviews [9]. 

Sonoluminescence is a complex phenomenon, involving as it does the motion of the 
bubble, the dynamics of the gas inside the bubble, and the mechanism that produces the 



flash of light. Our main concern in this work wiU be the last of these. The literature 
contains two classes of models to explain the flash of light. One involves the gas inside the 
bubble in an essential way [10]. There is no doubt that the gas undergoes compression and 
heating during the contraction phase of the bubble's motion, and this type of explanation 
relies on either thermal radiation, or else bremmstrahlung, to produce the light. 

The second type of explanation, on which we shall focus, is that the observed light 
is due to "vacuum radiation" [11-15], which is a dynamical counterpart to the well-known 
Casimir effect. In this view, given a particular motion of the bubble, 

r = R{t) (1) 

(here r is the radial coordinate that describes the bubble's surface, and R{t) is a prescribed 
function of time; we assume a spherical bubble centered at the origin for simplicity), the 
radiation would take place even if the bubble were completely evacuated. The role of the 
gas, and in particular the special role that seems to be played by the noble gases, is merely 
to modulate the motion of the bubble, i.e. to give rise to a specific R{t). It is then the 
motion of the boundary that directly gives rise to the radiation. 

To explore this idea, our approach will be to take as given all of the physics associated 
with the motion of the bubble and the dynamics of the gas, and to extract therefrom the 
single function R(t) which represents the experimentally measured bubble motion. Our 
next task is to construct a model in which the electromagnetic field is coupled to the bubble 
surface at r = R{t). 

One approach would be to attempt to derive this coupling from an examination of 
the dielectric properties of water; in practice, this would mean simply endowing the water 
with a dielectric constant e, and letting 

e(f , t) = ee(r - R(t)) + 9{R{t) - r). 

Considerable attention in the literature has been devoted to the case of R{t) — const., 
which, although it obviously neglects the dynamical mechanism that turns Casimir energy 
into real photons, is supposed to provide an order of magnitude estimate of the energy 
available, which can then be compared to the energy that is produced in sonoluminescence. 
The question of whether the Casimir energy is sufficient in this respect has become a rather 
controversial one [18]. Attempts to treat this problem dynamically have led to interesting 
results, but have not fully resolved the issue [12,13,14,19]. 



In this work, we shall choose a coupling that is not derivable (at least by us) from 
a direct consideration of the underlying physics. Rather, the interaction is chosen both 
for its simplicity and because it naturally leads to a coupling localized on the boundary 
r = R{t). In addition, as well shall show below, when one considers the case of two static 
parallel plates (Casimir's original problem) one recovers precisely the original Casimir 
energy, and is therefore encouraged to hope that the model may be a valid representation 
of the dynamical situation as well. 

In section 2 we shall introduce the model and derive the boundary conditions on E 
and B that it implies. In section 3, we look at two instructive cases that are not directly 
related to sonoluminescence: the case of static, parallel plates mentioned above, and the 
case of a strictly time-dependent source, with no spatial dependence. In this latter case, 
we shall discover the existence of unstable modes that can lead to production of radiation 
at unexpectedly large rates. 

In section 4, we tackle the case of greatest interest, the collapsing bubble. Even 
classically, we are unable to solve the equations exactly, but we develop an approximation 
scheme that relies on the fact that the radius of the bubble is small, in the sense that 
R{t) « cT, where T is a time characteristic of the width of the sonoluminescent pulse. 
In this approximation we find the same sort of unstable modes that existed in the purely 
time-dependent case. Section 5 is devoted to conclusions, and we have also included an 
appendix in which further properties of the FF interaction term are discussed. 

II. The Model 

We consider the following Lagrange density [16]: 

£=-^[F^,F^^ + /(x)F^,Fn • (2) 

Here F^^, has its usual meaning: F^^, = d^A^y — d^A^ where A^ is the 4- vector electromag- 
netic potential, and F^jy = ^^^lypaF^'^ , where e^i^p^ is the totally antisymmetric symbol on 
4 indices, and e0i23 _ 

As is well-known, Ff^^F^^ is a total divergence: 

F^uF^'' = d^[2e^^P''A,dpA,] (3) 



so up to a surface term, C can be written 



(4) 



Hence, by choosing 

fix) = mgix)] (5) 

we obtain 9^/ = fod^g5{g) so the second term in C represents the coupling of the electro- 
magnetic field to the surface given by g{x) = 0. (/o is a dimensionless constant.) For the 
case of the sonoluminescing bubble, we would choose g{x,t) = r — R{t). 
The equations of motion that follow from C are 

d^[F^^'' + f{x)Fn=Q , (6) 

or, since = identically, 

^^F''^ + {^^f)F^^^ = . (7) 

— * — * 

If we define E and B in the usual way, we obtain the modified Maxwell equations 

V • F + V/ ■ 5 = (8) 

ixB-E-fB-ifxE^O (9) 

— * 

together with y • 5 = and \/ x E + B = 0. With the choice (5), we have / = fQinP5{g) 

— * — * 

and v/ = fa^^id) where = {g, \/g) is the four-normal to the surface. To see what 
these equations entail, we write 



E = E^e{g) + E2e{-g) (10a) 
B = B^d{g) + B^ei-g) (106) 



and substituting into equations (8) - (9) , we find that the pair {Ei,Bi) satisfy the free 
Maxwell equations for ^ > 0, and likewise {E2, B2) satisfy them for ^ < 0. At = 0, we 
have the boundary conditions: 



n-{Ei- E2) +fon-B = 



(11) 



n X (Si - B2) - no{Ei - E2) - fo{noB + n x E) = 



(12) 



n-{Bi- B2) = 



(13) 



fix {El- E2) + no(5i - S2) = . 



(14) 



Notice that the second pair of equations, (13) - (14), removes the ambiguity as to which 



III. Special Cases 

(a) Parallel plates 

Before dealing with the time-dependent case, let us explore the physical significance 
of our model by revisiting the case originally considered by Casimir [17], i.e. two infinite 
parallel planes separated by a distance a. We take g{x) = z{z — a). The planes divide 
space into 3 regions, z<0, 0<z<a and z > a, and in each region we can choose a 
plane-wave solution to Maxwell's equations: 




(15) 



where k — {ki, /c2, ks) and k' — {ki, /c2, — fcs) and 



(16) 



We need both k and k' because the boundary conditions will mix them. 
Maxwell's equations imply that 



k xb = -oje k X e = ujb (17) 

and 

k' xb' = -ue' k' xe' =ujb' (18) 

in each of the three regions, which in turn imply the a;^ = = k'"^. 

It is now a matter of implementing the boundary conditions at 2; = and at 2; = a. 
After some algebra, it is not hard to show that the content of these conditions reduces to 

z = Q: 63 + = (19) 
z = a: b^e^^^" + b'^e-^^^"" = (20) 

(here 63 means the third component of b in the region < z < a, and similarly for 63) 
from which it follows that 

kza = mT, n = 0,±l,±2, . (21) 

This is exactly the same spectrum used by Casimir in his original paper, and therefore the 
Casimir energy 5E will be the same as his result: 

^ - (22) 

where L? is the area of one of the plates. 

(h) Time dependent source 

Armed with the knowledge that our model reproduces the static Casimir energy, we 
now proceed to another simple example, which is very different physically: we take /(x, t) 
to depend only on t: 

/(x,t) = 0, t<0 

= gt, 0<t<T (23) 
= gT, t>T . 



Since / is not a step function, there is no bubble in this case. Rather, the source is turned 
on everywhere at once at t = 0, and is turned off again at t = T (/ = const, is without 
physical consequence, because F^^^F^^'^ is a total divergence). 

The equations of motion in this case are of course just the free Maxwell equations for 
t < and t > T, whereas for < t < T, one of the Maxwell equations is modified: 



^ xB-E^gB . 



(24) 



One can now study a plane-wave solution, propagating, say, along the ^-axis. For t < 
we write 



B = {ax + 6y)e^('=^-'"*> 
E = {bx - ay)e'^'''-'^'\ 



(25) 



with k'^ = uj'^. At t = 0, this will be matched to a solution of the form 



B = e 



ikz 



e-'^\ax + (3y) + e*"*(7£ + 5y) 



iClti 



E = e 



ikz 



e-'^\(3x - ay) + e'''\5x - ^y) 



(26) 



where, because of the modification to Maxwell's equations, one has 0,^ = fij. = [k{k± g)]. 
Corresponding to each of these solutions is a particular polarization C± — -^[xTw]- 

— * 

When matched to the t < solution, the expression for B, < t < T, becomes 



20 



+ 



V2 



20 



+ 



V2 



.k + 0_ a — ib 

+ {—rr: ) —C^e 



20_ 



V2 



20, 



V2 



(27) 



with a similar expression for E. 

One can extend this analysis by matching this solution to a suitable expression for E 
and B in the region t > T, where of course a;^ = /c^ again. But we shall not need this 
extension in what follows. 

The feature most worthy of note is that, (for ^ > 0) the frequency 0_ becomes 
imaginary when k < g. (If ^ < 0, then 0+ becomes imaginary.) Hence E and B grow 



exponentially with time over the interval <t <T. We shall see below that, at least in a 
certain approximation, this feature persists in the case of a spherically oscillating bubble. 

— * — * — * 

To quantize this model, we can express E and B in terms of a vector potential A, and 

— * 

endow the fourier coefficients of A with the appropriate commutation relations. Effectively 
this means that the coefficients a and b in the above expressions become quantum operators. 
We must also generalize our solution to the case of a plane wave propagating in an arbitrary 
direction, but this is easily done since the z-axis used above was in no way special. It is of 
interest to compute the rate of energy production per unit volume by the external source. 
We do this by forming the hamiltonian density 

n=l{E-E + B-B) (28) 

and taking the vacuum expectation value of its time derivative. 'H is normal-ordered so 
that (0 I I 0) = for t < 0. Since the expressions are quite lengthy, we simplify matters 
by retaining only those pieces which grow exponentially. After some calculation, we then 
find 

^(0 I « I 0).., = mOiT-i)^,ll e^^^V^U. (29) 

where the notation "exp" on the matrix element means the exponentially growing piece. 
A simple stationary-phase estimate of the integral gives 

j^{0\n I 0)e.p e{t)9{T -t)-£^ e^' . (30) 

We can try to connect this to sonoluminescence (despite the fact that there is no bubble) 
by choosing T = lO"-*^^ sec (the duration of a typical pulse) and 1/g = 2x 10~^ m. (the 
cutoff on the observed spectrum). We then find gT = 1.6 x lO'^, which, needless to say, 
produces a huge number when inserted in the exponent in eq. (30). 

At this point, we can simply argue that our model is too far removed from the phe- 
nomenology of sonoluminscence to be expected to give reasonable results. Later, however, 
we shall have to deal with this question in the context of a more realistic model, to which 
we now turn. 



IV. The collapsing bubble 

We take f{x) = fo9{r — R{t)). Our strategy will be to attempt to solve the classical 
problem, looking for the kind of exponential behavior in time that we found in the previous 
example. If this is indeed found, then, reasoning by analogy with the previous example, 
we will argue that, when quantized, the model will produce an exponentially growing pulse 
of vacuum radiation over some period of time. 

Because of the spherical symmetry, it is appropriate to expand E and B in terms of 
vector spherical harmonics. Different values of / and m will not couple to each other. 

We write 

E = eiLYi^ + (re^ + e2)v>lm + ^^^^earK/m (31) 

and 

B = biLYim + irb'^ + b2)vyim + ^^^^^2^=^/^ . (32) 

Here L = jfx\/, and the e's and 6's are functions of r and t; we cannot separate variables 
any further because the boundary conditions will mix r and t. 

These forms automatically satisfy \/ ■ E = sy ■ B = 0. The rest of Maxwell's equations 
imply that 

0-^"-l^'+^l±ll^ = O. (33) 

where (p is any of the e's or 6's, and furthermore that 

ei = ib^ and bi = —ie2 ■ (34) 

In these equations, a prime is ^ and an overdot is ^ . So if we cast the boundary conditions 
entirely in terms of 62 and 62, we can recover ei and bi from equation (34). 

In fact it is not hard to express the boundary conditions in terms of 62 and 62. We 
expand E and B as in eqs. (31) and (32) separately inside and outside the bubble, and we 
let A(fi = (fiout — <fin where once again (f is any of the e's or 6's. Then we find, at r = R{t): 



A62 = ; 



(35) 



Ae2 = -fob2 



(36) 



A62 = [| + ^2 + Re2] (37) 

and 

Ae2 = + ^2 + ^^2] • (38) 

One can show that the expressions on the right-hand side of these equations all have zero 
discontinuity at r = R{t), so there is no ambiguity as to which values to insert. 

For simplicity, we choose to analyze the case 1 = 1. Then the most general solutions 
to equation (33), for 62 and 62, are: 



and 
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(39) 



Here the /9's and 7's arc arbitrary functions of the indicated arguments. In writing these 
equations, we have imposed the requirement that 62 and 62 be regular at r = 0. The 
functions /3 and 7 determine the waves propagating inward from infinity and should be 
taken as initial data. In principle, it should be possible to use the boundary conditions, 
eqns. (35) - (38), to determine the inside solution, specified by 13 and 7, and the outgoing 
waves specified by 13 and 7. The effect we are looking for is to see whether non-exponentially 
growing incoming data (3 and 7 can generate exponential growth in the outgoing solution 
/3 and 7. 

Unfortunately, when we substitute the forms (39) into the boundary conditions, we 
find rather complicated functional difference equations that we do not know how to solve. 
Instead we rely on an approximation that is based on the following observation. Phe- 
nomenologically, the smallest time scale we are interested in is the width of the sonolu- 
minescent pulse, At^j„ = lO"-^-*^ sec. The largest length scale we are interested in is the 
maximum size of the bubble, Rmax — 10~^ m. But in units with c = 1, 



At. 



max 



rmn 



10 m/sec 



30 



1 



(40) 



(this is actually a generous overestimate, since sonoluminescence occurs when the bubble 
radius is at least an order of magnitude smaller than Rmax)- Thus it might make sense to 
regard R{t) as a small parameter, and to expand our equations accordingly. We do this as 
follows: in the expressions for 62 and 62, equation (39), we replace the arguments t±r by 
t ± er, where e is a bookkeeping parameter in which we perform a systematic expansion. 
At the end we set e = 1. 

To obtain a consistent expansion, we not only expand the arguments of the functions, 
we must also expand the functions themselves: 



where (p stands for any of the unknowns, 7, (3 and 7. The input functions (3 and 7 are 
regarded as known and are not so expanded. The advantage of this expansion procedure 
is that we obtain thereby relations among functions all of which are evaluated at the same 
argument t. Henceforth we denote ^ by a prime. 

To obtain non-trivial results, we must retain terms up to order e^. For convenience, 
we introduce the notation pi{t) = P'o'it) and P2{t) = 7o'(*)- We find 



(p{t) = (po{t) + eipi{t) + e^<p2{t) + • • • 



(41) 



Po = P ; 70 = 7 



A = /32 = 7i = 72 = . 



73 = ^i?' [P2 - 1"' - foPl] . 



(42) 



and 




) [(1 - 3i?2) + RRp[] . 



) [2p2 + RRp'2] 



(43) 



The functions whose behavior we want to study are $3 and 73, which give the first non- 
trivial corrections to the outgoing waves $ and 7. Our strategy will be to solve eqn. 



(43) for pi and and then evaluate $s and 73 from eqn. (42). If we make the further 
approximation « 1, which must surely be true for any realistic bubble motion, we 
can drop the terms on the r-hs of eqn. (43), which can then be rewritten as 



- foRR fp[ 



(44) 



where M is the matrix { ^"^^ -2^r 1 • The eigenvalues of M are 

-1 ^/o ' ^ 



1 



1 



(45) 



Note that A_|_ > 0, A_ < 0. Let 9 be the orthogonal matrix 



V2 



/c_|_ —k- 
kj^ 



(46) 



where A;± = [1 ± Then, setting ^ ~ ^ (^p^ ) ' equation becomes 



-foRR 



^1 



X- (72 



+ 



K2 



(47) 



where 



«2 



= 6* 



/// 



Thus we have to solve 



-foRR I 

a = Act + K 



(48) 



whose solution is 



ait) = exp[- ! dt ' -^] - I dt' exp{ f 
J to foRR . J to J to 



3A 3k 



dt"- 

'to foRR foRR 



+ crito) 



(49) 



There are two such solutions, one for A+ and one for A_. Generally, the exponential 
behavior that is manifest on the r.h.s. of eqn. (49) will cancel in the first term, but will 
survive in the term proportional to cr{to). Because in the region of interest we have RR < 0, 
it will be A+ that gives the exponentially growing behavior (for /o > 0). 



How do we interpret this result? First, we must recognize that the approximation we 
have made is potentially very dangerous, because the highest derivative in eqn. (48) is 
multiplied by a factor, ^ which we expect to be quite small, and indeed which we 

expect to go to zero for | t | large. [Here we are restricting ourselves to only one cycle of 
the bubble's motion, so we take R{t) — > const, as | t |— * oo.] As a consequence, it appears 
from the solution that the exponential behavior becomes more pronounced the smaller RR 
becomes, whereas we know from the original equation, (48), that for RR strictly zero the 
solution is just a = —jf^, which exhibits no exponential behavior at all. 

Our response to this is to imagine that for | foRR \ below some threshold value, it is 
indeed negligible, and therefore a — —jK. At some time to, \ foRR \ crosses the threshold, 
and the solution (49) kicks in. We therefore fix the arbitrary constant (j{to) to be —jK{to). 
As we have already observed, for one of the two choices of A, a will be an exponentially 
growing function (except for the unlikely possibility that K,{to) = 0). 

The good news is that, within the context of our approximation, we have found the 
exponential behavior that we are looking for. The bad news is that, just as in the earlier, 
simpler example, we have an embarrassment of riches. There appears to be no mechanism 
within the model for turning the exponential behavior off - we have already fixed the one 
free parameter cr{to). Let us assume that there is a dissipative mechanism, having to do 
with the properties of the gas inside the bubble or the liquid outside it, that 

we should add to our model. This will abort the vacuum radiation after some charac- 
teristic time. Because the exponential growth produces so much radiation, we must assume 
that the abortion takes place after only about a femtosecond. In this picture, therefore, 
the observed pulse really consists of an exceedingly rapid exponential rise of duration a 
femtosecond or so, followed by a relatively slow fall lasting tens of picoseconds, for which 
the dissipative mechanism is responsible. 

V. Conclusions 

It has been argued in the literature that vacuum radiation cannot be the source of 
sonoluminescence because the static Casimir energy is so small. Without entering into 
the controversy over how big the static Casimir energy is for a spherical bubble [18], we 
believe that, as illustrated by our model (which, after all, correctly reproduces the static 
Casimir energy in the case of parallel plates) there are two additional factors that ought 
to be taken into account: 1) the Casimir effect arises essentially from the coupling of the 
electromagnetic field to a boundary. When that boundary is moving, the field is coupled 



to a time-dependent source, which in and of itself leads to the production of energy; and 
2) if this time-dependent coupling gives rise to unstable modes, as it does in our model, 
then an unexpectedly large amount of energy can be produced. 

The present work raises a number of issues for further investigation. Perhaps most 
important is tightening up the approximate treatment we have given for the classical 
solutions in the case of the collapsing bubble. It would be preferable to have a method 
of analysis that would conclusively demonstrate whether the period of exponential growth 
exists, and whether the model contains not only a mechanism for turning on the pulse 
but also for turning it off. Failing that, it will probably be necessary to include additional 
physics having to do with the kind of dissipation mechanism discussed above, capable of 
damping the vacuum radiation to a level compatible with what is seen experimentally. If 
this is the case, then the shape of the sonoluminescent pulse should exhibit a very rapid 
rise followed by a much slower decay. 

Within the context of the expansion employed in this paper, as might be expected 
one finds that the higher orders in e become progressively more complicated. We have 
examined the next non-trivial term, which is e^, and we have verified that it does not 
qualitatively change the exponential behavior found in order e'^. We have not checked, 
however, that the contribution is numerically small compared to the contribution. 

It would be useful as well to quantize the electromagnetic field in the presence of 
the collapsing bubble, much as we did for the simpler purely time-dependent case. We 
believe that any exponential behavior in the classical system will persist in its quantum 
counterpart, but having the explicit expression for quantum vacuum radiation would allow 
one to compare the details of the photon spectrum with experiment. It would also be 
useful, for numerical work, to have a good analytical approximation to R{t). 

Ultimately, electromagnetic radiation can be produced only by charges in motion. In 
the case of sonoluminescence, whether those charges effectively reside at the boundary of 
the bubble, as we contend in this paper, or within the gas inside the bubble, is a question 
that still awaits definitive resolution. 
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Appendix: More about FF 

In this Appendix we add a few remarks about the properties of FF. The f{x)FF 
interaction analyzed in the text was chosen as the most convenient form for couphng the 
electromagnetic field to the bubble boundary. Whether it correctly captures the essential 
physics of this coupling is a matter that will require further investigation. 

Another way of writing FF is just 2E ■ B. It is, apart from the familiar F^i.F'^'^ = 
— ^^), the only Lorentz invariant that can be constructed from E and B by algebraic 
means. A term in the Lagrangian of the form OF^^F^'^ , with 6 constant, has no physical 
consequence in an Abelian gauge theory such as electrodynamics, because, as noted in eq. 
(3), it is a total divergence. In a non- Abelian theory (such as quantum chromodynamics 
or the electroweak theory) this term is still a divergence, but it nevertheless gives rise to 
non-perturbative physical effects because of the non-trivial topological structures, called 
instantons, that exist in such theories. 

The anomalous divergence of the U{1) axial vector current is proportional to FF. 
This term is directly responsible for the decay of the 7r° meson into 2 gamma rays, which 
is its dominant decay mode. 

In string theory, because of a property called ^'-duality [20,21], two seemingly different 
theories can in fact be equivalent. This can be very useful because often one of the two the- 
ories is strongly coupled (and therefore interactable) whereas the other is weakly coupled. 
As it turns out, the system we have been studying in connection with sonoluminescence 
exhibits a simple form of S'-duality. To see this, it is useful to consider a slightly more 
general Lagrangian: 

C = -i[^i(x)F^,F^^ + ^2(a:)F^.F^^] . (Al) 

In the text we had ipi = 1 and <^2 = f{x)- In string theory, ip\{x) would be related to the 
dilaton, whereas (p2{x) is known as the axion field. 
This system is governed by two sets of equations: 

d^F^'' = (A2) 

and 



(^3) 



If we express F^j^ in the usual way as 

F^. = d^A, - d^A^ , iA4) 

then (A2) is an identity, while (A3) is a dynamical equation obtained by varying the 
Lagrangian with respect to Af^. 
We define 

G^i/ = ^^iF/^i^ + V^Ff^i. ■ {A5) 

In Minkowski space, as is easily shown, the dual of a dual is the negative of the original 
tensor. Therefore 

G^u = -^iFfj,„ + (p2Ffj,,y . (AG) 

We can invert these relationships to obtain F and F in terms of G and G, and then 
substitute them into C: 

^ = -^(^^) {-fiG^uG^'' + ^iG^.G^^\ . (A7) 

If we set 

and vary C with respect to S^, we obtain 
and also 

d^G^"" = (AlO) 



as an identity. It is straightforward algebra to show that (A9) is the same as (A2), and 
(AlO) is, by definition, the same as (A3). Thus the physical content of the Lagrangian (A7) 
is the same as (Al), but the dynamical equation in one case is an identity in the other case. 



and vice versa. We see that if we can solve a system with sources (991, (^2), then by duahty 
we automaticaUy obtain a solution with sources (—2^^^, -^^r-^)- It is not clear, however, 
whether practical use can be made of this observation in the case of sonoluminescence. 
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